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The elegance and usefulness of a complex formulation of the basic lensing equations is demon- 
strated with a number of applications. Using standard tools of complex function theory, we present, 
for instance, a new proof of the fact that the number of images produced by a regular lens is always 
odd, provided that the source is not located on a caustic. Several differential and integral relations 
| between the mean curvature and the (reduced) shear are also derived. These emerge almost au- 

■ tomatically from complex differentiations of the differential of the lens map, together with Stokes' 
^\ ' theorem for complex valued 1-forms. 
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Gravitational lensing has become one of the most important fields in present day astronomy. The enormous activity 
in this area has largely been driven by considerable improvements of observational capabilities. Gravitational lensing 
' has the distinguished feature of being independent of the nature and the physical state of the deflecting mass. It is 
[ therefore perfectly suited to study the dark matter in the Universe 0,0. 
t-H ■ One of the issues which has recently attracted a lot of attention is concerned with parameter-free reconstructions of 
projected mass distributions from weak lensing data. (For a recent review, see j3j.) Thanks to new wide-field cameras 
and imaging with 8?7j-class telescopes, the quality of the data is expected to increase rapidly. Initiated by a paper 
of Kaiser and Squires B , a considerable amount of theoretical work on various reconstruction methods has recently 
also been carried out H], ^ . The main problem consists in the task to make optimal use of limited noisy data in a 
' parameter-free manner, that is, without modeling the lens. 

In the present paper we take up some of the theoretical discussions and demonstrate rather systematically that the 
Q ■ complex formulation of lensing theory often simplifies things considerably. In particular, a number of equations which 
^—j | are used in mass reconstructions, emerge almost automatically. 

In outline, the paper is organized as follows: For reasons of self-consistency, we provide in Section 2 a brief derivation 
of the basic lensing equations that are used in the remainder of the paper. These are then translated in Section 3 into a 
complex formulation, where some mathematical tools are recapitulated as well. It will turn out that the reconstruction 
problem is basically equivalent to the task of solving the so-called Beltrami equation, at least for noncritical lenses. 
This part of the paper has considerable overlap with JtJ. Turning to applications in Section 4, we give - as far as we 
know - a new proof of the fact that for a regular lens the number of images is always odd, provided that the source 
is not located on a caustic. The proof uses only standard tools of complex analysis, which are, for instance, familiar 
from derivations of the theorem of residues. One of these formulas is an explicit expression for the index of a closed 
curve relative to a given point. Next, we derive several relations between the mean convergence and the (reduced) 
shear by (repeated) applications of the complex differential operators d/dz and d/dz to the differential of the lens 
map. Several other useful relations for lensing reconstructions, involving integrals over bounded domains, are derived 
at the end of the paper. 

The purpose of this article is mainly methodological. We hope that others will take advantage of it, especially in 
teaching the pleasant field of gravitational lensing. 



II. BASIC LENSING EQUATIONS 

For the benefit of those readers who have not studied the extensive monograph of Schneider, Ehlers and Falco [|l|, 
we start by giving a brief derivation of the basic lensing equations. 

The conceptual basis of gravitational lensing theory is extremely simple. This is at the same time one of the main 
reasons why it is so important for the astronomical study of mass distributions on all scales. For all practical purposes 
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the ray approximation for light propagation is sufficiently exact. In this limit the rays correspond to null geodesies in 
a given gravitational field g, and the evolution of the polarization vector is governed by the law of parallel transport. 
(These laws can be deduced from Maxwell's equations ||.) The null rays are orthogonal to the surfaces of constant 
phase, {S = const}, where S is subject to the eikonal equation 

g^d^SduS = 0. (1) 

For sufficiently strong lenses the wave fronts develop edges and self-intersections. Clearly, an observer behind such 
folded fronts sees more than a single image. This is the region of what is called strong lensing and occurs astronomically 
only rarely. 

Here we restrict ourselves to almost Newtonian, asymptotically flat situations. Generalizations to the cosmological 
context are easy and basically amount to interpret all distances in the formulas given below as angular distances. (For 
details we refer again to 0, hereafter quoted as SEF). The metric is then given by 

g = (1 + 2C7) dt 2 - (1-2C) dx 2 , (2) 

where U is the Newtonian potential. The spatial part of a light ray satisfies Fermat's principle, 

-0, (3) 



/goo 

for variations with fixed end points || . Here da 2 denotes the spatial part of the metric (||) . 

All this can be summarized by saying that gravitational lensing theory is just usual ray optics with the refraction 
index 

n(x) = 1 - 2U(x) . (4) 



In particular, the ray equation holds, 



dx 
ds V ds 



= Vn, (5) 



where s is the euclidean path length parameter. (Since light deflection is a scattering process, we can from now on 
forget about non-euclidean geometry.) 

In terms of the unit tangent vector e = dx/ds, eq. (^) can be written in sufficient approximation as 

— e = -2V x U, (6) 

as 

where Vj_ denotes the transverse derivative, Vj_ = V — (e, V)e. This gives for the deflection angle a — ei n — e/» n , 
with initial and final directions ej n and efi n , respectively, 

V ± Uds, (7) 

where the integral is taken over the unperturbed path (u.p.). Here, we insert the expression for the Newtonian 
potential of a mass density p(x). In the well-justified approximation where the extension of the lens (for instance a 
cluster of galaxies) is much smaller than the distances of the observer and the source to the lens, one finds readily 
(SEF, Chapter 4) 



««) = 4 G L, rr s (0 d 2 t' - ( 8 ) 



M 2 I £ - £' I 2 

where £(£) denotes the projected mass density on the lens plane. (For a point mass this reduces to Einstein's 
prediction of light deflection.) 

Combining this with elementary geometry, we arrive at the lens map for a given From Fig.l, which summarizes 

the notation of SEF, we read off the lens equation 

ri = ^£ - £>*&(€), (9) 
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which defines a map from the lens plane to the source plane. 

It is convenient to write this in dimensionless form. Let £o be a length parameter in the lens plane (whose choice will 
depend on the specific problem), and let r]o be the corresponding scaled length in the source plane, 770 = (D s /D d )£,o- 
We set x = £/£ci) V = ri/vo an d (following SEF) 



£(£oa0 , s D d D ds A 
K ( x ) = i "W = c n oc{£ x), 



with 



(10) 



^ r.r it — 
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Then eq. (m reads as follows 



y = x — ot(x) , 



(11) 



whereby eq. (H) translates to 



a(x) 
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— K{X>) . 



IT J]R 2 I X — X' 

It is obvious that a is a gradient of a two-dimensional Newtonian potential: 



(12) 




A A 



D. 



Observer 



FIG. 1. Notation adopted for the describtion of the lens geometry. 
a = V?A , ^ = 2G * k , G(a;) = — In | x | (* denotes convolution) . 



(13) 



Since G is a fundamental solution of the two-dimensional Laplace operator, ip satisfies the two-dimensional Poisson 
equation 



2k. 



(14) 
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For the differential Dip of the map p : 1R 2 — > 1R 2 , defined by eq. (|TT|), we use the standard parametrization 

D<p = f 1 -^^ "72 \ (15) 

V -72 1 - « + 71 / 

in terms of the mean (Ricci-) curvature k, determined by the trace of Dip, and the (Weyl-) shear vector 7 = (71,72). 
The eigenvalues of the symmetric matrix Dp are 1 — K=p | 7 |. The critical curves, satisfying det (Dp) — 0, are given 
by 

(l-nf- | 7 | 2 = 0. (16) 
The caustics are the images of these critical curves. In the vicinity of a caustic the amplification jj,, given by 

M = I dot (Dp) I ' 

becomes very large. 

In passing, we note that the lens map ( pi] ) can also be written as 

V a <f> = 0, with <j>(x,y) = \(x-yf - #c) • (18) 
This reflects the Fermat principle. Indeed, the delay of arrival times is directly given by the Fermat potential (f>: 

At = ^irFr^v)' ( 19 ) 

Examples of various lens maps are discussed extensively in Chapter 8 of SEF. Two standard cases are (with suitable 
choices of £ ) : 

Schwarzschild lens: y = x — x/\ x \ 2 ; (20) 
singular isothermal lens: y = x— x/\x\. (21) 

It is worth recalling the following general fact: In 1955, in a pioneering work of modern singularity theory, H. 
Whitney || studied generic properties of smooth mappings of the plane into itself and proved that the subset of 
mappings which have only fold and cusp singularities contains an open and dense set (with respect to the Whitney 
topology). Moreover, those maps of this set which satisfy a few mild global conditions are also stable. Clearly, these 
results are highly relevant to gravitational lensing. For realistic lenses we will only have folds and cusps, and no 
singularities of higher order. 

III. COMPLEX FORMULATION 

In this section we translate the basic lensing equations into a complex formulation. It will turn out that this is not 
only elegant, but also quite useful, because one can then apply various tools and techniques of complex analysis. This 
has also been noted before by other authors [Q. 

A. Mathematical Preliminaries 

We use standard notation when identifying 1R 2 with (F, by writing z = x + iy for (x, y) £ 1R 2 and dz — dx + idy, 
dz — dx — idy for the corresponding basis of 1-forms. In terms of the Wirtinger derivatives, 

a= d_lfd_.d\ d_ _ 1 / ' d_ .d_\ 

z = d~z ~ 2\dx~ ~ l dy~) ' z ~ dz ~ 2 [dx + % dy) ' ( } 

the differential of any smooth complex function / on (E has the representation 

df -fz dZ + % d ~ Z - ^) 



4 



We shall also write f z and / 2 for d z f and dgf, respectively. A function / is holomorphic if and only if dgf = 0. In 
terms of the Wirtinger derivatives the Laplacian is given by 

A = 43,0*. (24) 

We shall make repeated use of Stokes' theorem for complex- valued differential forms ond (or an open subset): If 
£1 is a compact subset of (F with smooth boundary dQ, then for every complex differential 1-form u> 



du = u. (25) 
n Jen 

An immediate corollary of eq. (|2^) is the Cauchy-Green formula: For a smooth function / we consider 

" = /-^7. (26) 
and apply Stokes' theorem (]2jJ) for fl minus an £-disk with center £. In the limit e — ► we obtain 

/(C) = ^ / + ^ / ^dzAdz. (27) 

For holomorphic functions the second integral is absent. (Note that dz Adz — ~2idx A dy.) 
The dilatation or Beltrami coefficient v — Vj of a smooth function / is defined by 

h ="ffz, (28) 
and this equation is also called Beltrami equation. Since the Jacobian Jf of / is given by 

Jf H/,1 2 - \h I 2 , (29) 

we conclude that \vt |< 1 if / preserves orientation and z// = if and only if / is conformal. For the interpretation of 
Vf we consider the infinitesimal ellipse field by assigning to each zed 1 the ellipse that is mapped to a circle by /. As 
indicated in Fig. 2, the argument of the major axis of this infinitesimal ellipse is [n + arg(uf)] /2, and the eccentricity 
e is 

e -\M + \h\-i + Wf\- m 




FIG. 2. Geometrical interpretation of the Beltrami parameter. 



Solving the Beltrami equation fl28| ) is then equivalent to finding a function / whose associated ellipse field coincides 
with a prescribed v. We shall see that this is just the inversion problem in gravitational lensing. Weak gravitational 
lensing corresponds to quasiconformal maps. A smooth map / is k-conformal if its Beltrami parameter Uf satisfies 
\i/f\ < k < 1. This means geometrically that there is a fixed bound on the stretching of / in any given direction 
compared to any other direction. 
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We now quote an existence and uniqueness theorem for the Beltrami equation. For a fixed k with < k < 1 let 
L oc (k, R) denote the measurable functions on(F bounded by k and supported in {z G <T\ \z\ < R}. 

Theorem : For v G L°°{z, R), there is a complex function f on(E , normalized so that f{z) = z + 0(1/ z) at oo, with 
distributional derivatives satisfying the Beltrami equation f% — v f z , and such that f s and f z — 1 belong to LP for a 
p > 2 sufficiently close to 2. Any such f is unique. The solution f is a homeomorphism of(U, which is holomorphic 
on any open set on which v — 0. If v G C 1 and v z G C , then f G C 1 . 

A proof of this theorem can, for instance, be found in JTo{ ] . 

The reconstruction problem (for noncritical lensing) will lead to the inhomogeneous Cauchy-Riemann equation 

d-J = h. (31) 
In case the smooth function h has compact support, the Cauchy-Green formula ( |27| ) provides one solution: 

f(0 = ^~ (J^-dzAd-z. (32) 

Obviously, / is only determined up to an additive holomorphic function. If the solution is assumed to be bounded, / 
is unique up to an additive constant. 

From the solution ( |32| ) we see that (ttz)^ 1 is a fundamental solution of the differential operator <9 2 , 

(33) 

because (|32] ) can be written as 

f = --*h. (34) 

TT Z 

A special case of the so-called Dolbaut Lemma in several complex variables implies that one may drop the assumption 
that h has compact support: 

Theorem : For any smooth function h ond there exists a smooth function f such that fts\ ) holds. 
For a complete proof, see Chapter 2 of 
As an easy consequence we have the 

Corollary: For any smooth function h there exists a smooth solution of the Poisson equation A/ = h. 
In the following we often use the abreviations d = d z , 8 = d z - 

B. The complex Lens Mapping and its Differential 

The lens mapping ip : M 2 i — ► M 2 , 

y = <p (x) = x — W (x) , (35) 
is now written as / : <T i — > (F, w = f{z) with z = x\ + ix2, w = y\ + iy^. We have 

f(z) = z- 2dip (36) 

or 

f = B(zz- 2i/>) . (37) 

Eq. ([l4] ) becomes 

2ddip = k . (38) 
The differential of / will be very important. From (^) and ( p8| ) we obtain 

df = (1 - K )dz - 2d 2 i)dz. 



G 



But 

0V = \ (?l - dl) V> + = \ (71 + *72) , 

according to the original definition (JlJ) of the shear vector. Introducing the complex shear 

7 = 7i + «72 (39) 

we obtain 

df = (1 - n)dz - 1 dz . (40) 
Hence, the Beltrami parameter Vf of the lens map is given by 

»f = ih- (41) 

This agrees with the reduced shear introduced by Schneider and Seitz [Q. 
The examples ( p0| ) and (^T|) become: 

Schwarzschild lens: f(z) = z — — , //<■ = — — ; (42) 
singular isothermal lens: /(z) = z — — , Vf — —. — r . (43) 



For reference, we note that the amplification [i is according to ( |17| ) , (|29|) and (^0|) given by 

= I J/l = II5/I 2 - \d]\ 2 \ = |(1 - nf - | 7 | 2 | . (44) 

IV. APPLICATIONS 

The usefulness of the complex formulation will be illustrated in this section with several applications. No new 
results are obtained, but some of the derivations become simpler and more natural. 

A. Number of Images for a regular Lens 

The important fact that the number of images for a regular lens is always odd, provided the source does not lie 
on a caustic, is traditionally proven with the help of some elements of Morse theory jjj. We now give a proof which 
uses only standard tools of complex function theory that are used, for example, in the derivation of the theorem of 
residues. In particular, we make use of the following analytic formula for the index of a closed (rectifiable) curve 7 
relative to a point 0^7: 

ind ( a ) = _L f . (45) 

Zwi J 1 z — a 

This index is equal to the winding number of 7 around a and hence an integer. Furthermore, it is a nomotopic 
invariant, changes sign under orientation reversion, and is additive under composition of closed curves (see, e.g., 
Chapter IV of @). 

Consider now a point w in the source plane with images f^ 1 (w a ) = {zi, . . . , zn} in the lens plane. The complex 
1-form 

— 1 ^ (46) 



2iri f — w 



is regular on (F\ \Jj D e (zj), where D e (a) denotes the closed disk with center a and radius e. It is also closed, and 
therefore Stokes' theorem (09) gives 
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If df A 1 f df 



2? ™ JdD R (0) f - w ° ~t 2?r * JdD t (zj) f -W' 

Now, for a closed curve 7 we have by the transformation formula of integrals and 

1 f df If dw 



27Ti J-y f — w 2iri 



— indfoy (w ) ■ (48) 

/07 — w o 



Asymptotically the lens map approaches the identity, and hence the left hand side of ( f47| ) is equal to 1 for R sufficiently 
large. Therefore, we have 

N 

1 = 22 indf odDe ( Zj -) (w ) = rii - n_i + 2 (n 2 - n_ 2 ) + ■ ■ ■ , (49) 

3=1 

where ua denotes the number of Zj in {zi, . . . , zjy} for which the index in (^9|) is equal to A. 

For the special case, when w Q is not on a caustic, the Jacobians Jf{zj) do not vanish and all indices are thus equal 
to ±1 (+1 if / is orientation preserving and —1 if it is orientation reversing at Zj). Hence 

N = n\ + ri-i , 1 = ni — n_i , (50) 

implying that 

N = 1 + 2n_i (51) 

is odd. 



B. Relations between mean Convergence and reduced Shear 

The Beltrami parameter (reduced shear) Vf of a lens map is in principle observable. What we are really interested 
in is, however, the mean curvature n which is related to the surface mass density by (|l0|). 
In view of ( p9|) it is natural to look first for relations between the complex shear 7 and k. 
Eq. ( ftp] ) for the differential of the complex lens map and (|3^) give 

7 = -df = 2d 2 4> . (52) 

In order to get a useful relation we differentiate (^||) and use (|3£ 

d 1 = 2d (dBtf) = 8k . (53) 

This can be regarded as an inhomogeneous Cauchy-Riemann equation for n. With the results in Section 3.1 we 
conclude 

7T \Z J 7T \Z J 

or 

1 1 fr.S 

n = ^*7 + k . (54) 

7T Z -2 

The additive constant n reflects the fact that a homogeneous mass sheet does not produce any shear ('mass sheet 
degeneracy'). The real form of (|54|) appears the first time in Q|. In making use of (pd]), we obtain an integral equation 
for k when v is known: 

k = - — \ * [v (1 - k)} + k . (55) 

7T Z z 

This has been used, for instance, in || for nonlinear cluster inversions. 

We add that (^4|) has an inverse, that also appeared in the influential paper p|] of Kaiser and Squires. From (|5 
and (|13|) we obtain 
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7 = 4<9 2 G*k. (56) 
Since the fundamental solution G of the two-dimensional Laplace operator is 

G=^\n\z\ = ^\n{zz) (57) 



we find 



7 = - — * « • (58) 



Note that ( p3[ ) has the real form (k is real) 



Let us differentiate (p3|) once more 



giving 



( 9 l7 i + J 72 j (59) 
1 C172 + o 2 7i ' 



aa K = 9 2 7 , (60) 



Ak = Ad 2 [v(l - k)] , (61) 



from where we could again arrive at (|55j). The mass-sheet degeneracy is reflected in the following invariance property: 
Eq. (|6"l|), for given i/, remains invariant under the substitution 

k — > Ak + (1 - A) , (62) 

where A is a real constant E3 . 

We can use (53) in a different manner. First, we write this equation as 

3k = d \v(l — k)] = (1 — n)dv — vdn . 

This becomes simpler in terms of K := ln(l — n): 

dK - vdK = dv . (63) 

To this we add its complex conjugate. Noting that K is real, we obtain again an inhomogeneous Cauchy-Riemann 
equation, this time for K: 

dK = h(y) , (64) 

whereby the inhomogeneity 

h{v) = (1 - |z^| 2 ) 1 [dis + lJdU\ (65) 

is in principal observable. 

The real form of this equation was obtained by Kaiser |l5| and has often been used in the analysis of cluster data. 
The complex version appears also in . 

It should have become clear at this point that the complex formulation is also useful. The relations, derived in this 
subsection, emerge alsmost automatically by just applying d and 8 to the coefficients of the differential of the lens 
map. 



C. Other useful Reconstruction Equations 

Real lensing data are always confined to a finite field of the sky. Therefore, the solution of flS4| ) in the form (j32|), for 
example, involving an integration over all of(T, is in practice not very useful. One can, however, also obtain integral 
formulas in which only integrations over bounded domains occur. 

In order to arrive at these, we write the inhomogeneous Cauchy-Riemann equation in terms of differential forms: 
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d"g = uj . (66) 

Here w is a 1— form and we use the standard decomposition d = d! + d" of the exterior derivative, satisfying 

d! o d! = , d"orf" = 0, d'o d" + d" o d! = (67) 

(see, e.g., 0]). We make also use of the ^-operator, which is related to complex conjugations as follows: If a 1— form 
a is decomposed as a = ot\ + a2, where a± is of type (1, 0) and «2 of type (0, 1), then 



* a = i(a\ — Qfa) . (68) 



The following identities are useful: 



kk a = —a , ka = *a , 

<i* (ai + Q!2) = irf ai — id" a.2 , 

★d'<7 = id"g , -kd g — —id'g , 
d*dg = 2id'd"g = Agdx A dy , (69) 

where g is a function. 

Let now VL <Z(E be a bounded domain with smooth boundary dQ and A — |fi|. We show that <? minus its average g 
over fi, 

g=\ [ gdx A dy (70) 



.4 

can be represented in the following form 

g-g= j *n / ^. (71) 

The 1— form a in the integral is given by 
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in terms of the real Green's function H , defined by 

7 - 

A 



(72) 



Aff-i = -5, (73) 



together with the Neumann boundary condition on dfl. 

This is a consequence of Stokes' theorem. The integrand in (fn]) is 

*«Aw = *aA d"g = -d" (g * a) + 2fd" {kd" H) . 

By making use of ( |69| ) we obtain for the last term 

gd" (*d"H) = 2igd"d'H = -gAHdx A , 

while the first term is given by 

-d" (g ★ =-d(g* d"H) . 

Hence, 

■ka A u> = — / 5* + <7 — • 

This is just ( [n] ) since the last integral vanishes, due to the Neumann boundary condition for H . Formulas equivalent 
to (fn]) have been much used by S. Seitz and P. Schneider ||. 
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The starting point for the derivation of another useful relation is ( p8| ) in the form 

d (f — z) = —ndz — jdz , 
If we wedge this with dz and add the complex conjugate of the resulting equation we find 

ndz A dz = — d [k (zdz — zdz) + jzdz — ^zdz] . (74) 



Taking the average according to (|70| ) we arrive at 

('jzdz — jzdz) 
§ {zdz — zdz) 

where (•) denotes the average along the boundary dVL: 



(75) 



= §K (zdz -zdz) 
§ (zdz — zdz) 

For the special case of a disk D r we have along the boundary z = re lip , zdz — zdz = —2ir 2 dip, hence 

*=<«> + (7t> , (77) 
where 7 t denotes the tangential component of the shear 

It = 7i cos 2ip + 72 sin 2ip . (78) 
This relation is not new (see Ref. ||). Noting that 

1 [ r 

k= — 7t I n(r , (p)r dr dip , (79) 



and thus 

dR 



d\nr 

we can use (|77|) to obtain the interesting connection 



2 (k) - 2R , (80) 



= -2<7t>. (81) 



rflnr 

This has recently been used in an analysis of weak lensing data |^|. A useful integral form of it is, in obvious notation, 



/ r 2 \ f r2 dr 
R(r 1 )-R(r 1 <r<r 2 )=2 ( 1 1 1 / ( 7 t) — • 



(82) 



The left hand side of this equation is what Kaiser and Squires call the ^-statistics, C( r ii r 2)- One can use general 
weight functions for the average process || and try to optimize the choice for the detection of mass overdensities || . 
Note also, that the integral on the right in (|8^ ) can be written as 

1 

2^ 



r2 dr 1 f ( 1 \ 

(7*)- = -/ ^[-l)dxNdy. (83) 



We conclude by pointing out another appearance of a Beltrami parameter in lensing theory. An often used method 
for describing the shape of a galaxy image uses the second brightness moments 

Qi]=T7~ — [ I (x) (xi - Xi) (xj - Xj)d 2 x , (84) 
Norm J 

where / (x) is the surface brightness distribution and x is the center of light of the galaxy image. Regard now 
Q = (Qij) as a linear map of M 2 . If this is interpreted as a map z i — ► w (z) of (E it reads 
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w = ~ (On + Q 22 ) z + ~ - Q22 + 2iQi 2 ) z = l;trQ [z + X z] , (85) 

where 

(Q11 -Q 22 + 2zQ 12 ) . . 

x = t^Q • (86) 

X is called the complex ellipticity and is clearly just the Beltrami parameter of the map (p5|). The intrinsic brightness 
moments Qf) of the galaxy are defined corespondingly and it is easy to see that — Dip ■ Q ■ Dip, Dip being the 
differential (Jlq) of the lens map. The interpretation of x just given, allows us to find easily the corresponding relation 
between \ and x ' ■ One just has to compose the map ( pq ) on the right and on the left with the linearized lens map 

w = (1 — k) z — jz . (87) 

This gives readily 

( s ) _ -2v + x + v 2 X 



X 

with the inverse 



l + M 2 -29fi^x) ' 



2, + X W+^') 

X l + M 2 + 23fi(^( s )) ' [ ' 

A real derivation of these formulas is quite akward. They are used in applications by averaging over a set of galaxy 
images, together with statistical assumptions about the intrinsic ellipticity distribution (for instance (x^) = 0), to 
determine the reduced shear v of the lens map. Here, we just wanted to point out that X has the interpretation of a 
Beltrami parameter and that the relations (pq) and ( p9| ) are very easily obtained in the complex formalism. 
We hope that the reader will find other examples of such simplifications. 
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